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(3- decay properties of proton-rich odd-A and even-even Krypton isotopes are studied in the frame¬ 
work of a deformed selfconsistent Hartree-Fock calculation with density-dependent Skyrme forces, 
including pairing correlations between like nucleons in BCS approximation. Residual spin-isospin 
interactions are consistently included in the particle-hole and particle-particle channels and treated 
in Quasiparticle Random Phase Approximation. The similarities and differences in the treatment 
of even-even and odd-A nuclei are stressed. Comparison to available experimental information is 
done for Gamow-Teller strength distributions, summed strengths, and half-lives. The dependence of 
these observables on deformation is particularly emphasized in a search for signatures of the shape 
of the parent nucleus. 


I. INTRODUCTION 

One of the most exciting challenges in current Nuclear Physics is the understanding of nuclear systems under 
extreme conditions [Q. The opportunities offered by recent experimental work using beams of exotic nuclei and the 
corresponding theoretical efforts to describe them are of great interest, especially for nuclear structure physics and 
nuclear astrophysics (||. 

The decay properties and cross sections for nuclear reactions of radioactive nuclei are fundamental to understand 
various phases in the stellar evolution, including the energy generation, the nucleosynthesis, and the abundance of 
elements. Since this information cannot be determined experimentally for the extreme conditions of temperature 
and density that hold in the interior of the star, reliable theoretical calculations for these processes are mandatory. 
In particular, the decay properties of proton rich nuclei are fundamental to understand the rp process (rapid proton 
capture nucleosynthesis), characterized by very large proton capture reaction rates on proton rich nuclei ||. Of special 
importance in this context are the waiting points like ,2 Kr, where the rp process is inhibited and the reaction flow 
has to wait for the relatively slow /3-decay to continue. The total half-lives of the waiting points determine the speed 
of nucleosynthesis towards heavier nuclei as well as the isotopic abundances. 

Decay properties and nuclear structure are intimately related. It is clear that a precise and reliable description of the 
ground state of the parent nucleus and of the states populated in the daughter nucleus is necessary to obtain a good 
description of the decay, and vice versa, failures to describe the decay properties would indicate that an improvement 
of the theoretical formalism is needed. 

From a theoretical point of view the physics of exotic nuclei, characterized by very unusual ratios of neutrons 
and protons, can be considered as a test for the already well established models of nuclear structure that are used to 
describe stable systems. Since the parameters and interactions used in the usual shell model or mean field calculations 
are determined in order to reproduce the properties of known nuclei, they may not always be appropriate for use in 
the calculations of nuclei approaching the drip lines. 

Different microscopic models to describe the /3-strength are present in the literature. Models based on spherical 
single-particle wave functions and energies with pairing and residual interactions treated in Quasiparticle Random 
Phase Approximation (QRPA) were first studied in Ref. jj|]. Deformation was included in Ref. |fj, where a Nilsson 
potential was used to generate single-particle orbitals. Extensions including Woods-Saxon type potentials ||, residual 
interactions in the particle-particle channel (JtJ, Hartree-Fock (HF) mean field with separable residual interactions 
treated in Tamm Dancoff approximation ||, selfconsistent approaches in spherical neutron-rich nuclei |fj] and based 
on an energy-density functional 0 , can be also found in the literature. 

In a previous work 00 we studied ground state and /3-decay properties of even-even exotic nuclei on the basis 
of a deformed selfconsistent HF+BCS+QRPA calculation with density dependent effective interactions of Skyrme 
type. This is a well founded approach that has been very successful in the description of spherical and deformed 
nuclei within the valley of stability 0 - In this method once the parameters of the effective Skyrme interaction are 
determined, basically by fits to global properties in spherical nuclei over the nuclear chart, and the gap parameters 
of the usual pairing force and the coupling strength of the residual neutron-proton pairing force are specified, there 
are no free parameters left. Both the residual interaction and the mean field are consistently obtained from the same 
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two-body force. This is therefore a reliable method, suitable for extrapolations into the unstable regions approaching 
the drip lines. It is worth investigating whether these powerful tools designed to account for the properties of stable 
nuclei are still valid when approaching the drip lines. 

One possible way to establish the validity of the known approaches as well as the limits of their applicability, is 
a systematic investigation of nuclei covering the whole range from stability towards the drip lines. Exploration of 
series of isotopes moving away from the region of stability would fulfill these requirements allowing us to learn how 
the adequacy of the description evolves in progressively more unstable nuclei. 

we apply this formalism to the study of proton rich 


Following the same criteria as in our previous work [11 
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Krypton isotopes, including even-even as well as odd-A isotopes for the first time. The reasons why this is a region of 
special interest to study /3-decay have already been stressed in Refs. 00 - They include the large Q-values in proton 
rich nuclei |l!|, the competition of nuclear shapes j|,[l6|] that characterizes this mass region, and the possibility to 
approach systematically the N = Z isotope. Thus, we can test the validity of our formalism and look for discrepancies 
when approaching the drip lines. 

The paper is organized as follows. We first summarize in Sect. 2 the method of calculation. Mean field and QRPA 
with residual spin-isospin interactions are introduced and explained separately for both even-even and odd-A nuclei. 
Gamow-Teller strengths, excitation energies, Qec values and half-lives are also discussed. In Sect. 3 we present 
our results. We discuss similarities and differences in even-even and odd-A nuclei and analyze in detail the /3-decay 
observables in Kr-isotopes, comparing them with the available experimental data. Sect. 4 contains the conclusions 
and some final remarks. 


II. THEORETICAL FORMALISM 

In this Section we summarize briefly the theory involved in the microscopic calculations. More details can be found 
in Refs. 00- Our method consists in a selfconsistent formalism based on a deformed Hartree-Fock mean field 
obtained with a Skyrme interaction, including pairing correlations in the BCS approximation. We consider in this 
paper the force SG2 |l7j of Van Giai and Sagawa, that has been successfully tested against spin and isospin excitations 
in spherical 0]and deformed nuclei JI^]. Comparison to calculations obtained with other Skyrme forces have been 
made in Refs. |1 l| . [i~2| ], showing that the results do not differ in a significant way. The single particle energies, wave 
functions, and occupation probabilities are generated from this mean field. 

For the solution of the HF equations we follow the McMaster procedure that is based in the formalism developed in 
Ref. [jltl and described in Ref. |pc|| . Time reversal and axial symmetry are assumed. The single-particle wave functions 
are expanded in terms of the eigenstates of an axially symmetric harmonic oscillator in cylindrical coordinates. We 
use eleven major shells. The method also includes pairing between like nucleons in the BCS approximation with fixed 
gap parameters for protons A w , and neutrons A„, which are determined phenomenologically from the odd-even mass 
differences through a symmetric five term formula involving the experimental binding energies |2l| ]. The values used 
in this work are the same as those given in Ref. [Q. 

For odd-A nuclei, the fields corresponding to the different interactions were obtained by doing one iteration from 
the corresponding selfconsistent field of the closest even-even nucleus, selecting the orbital occupied by the odd 
nucleon according to the experimental spin and parity. For those cases where this experimental assignment is not well 
established we choose the orbital closer to the Fermi level. The chosen state is blocked from the BCS calculation, and 
we assign to it a pair occupation probability of 0.5. The effect of doing several more iterations from the even-even 
case, in order to see how the extra particle polarizes the core, was studied in Ref. [^2| without observing significant 
changes. We repeated these calculations here looking for some effect on the GT strength distributions, but again the 
changes were negligible. According to this procedure our spin and parity assignments are as follows: 5/2 + for 77 Kr 
and 75 Kr (as determined experimentally pl|), 3/2“ for 73 Kr from the most recent experimental determination p3| . 
For 69 Kr and 71 Kr we take the spin and parity according to our calculations and assign K v = l/2 _ in the first case 
and K n = 9/2 + or K 7T = 3/2“ in 71 Kr depending on the oblate or prolate shape. We describe the even Z odd N 
parent nucleus by removing one neutron from the selfconsistent field of the even-even nucleus and the odd Z even N 
daughter nucleus by removing one proton. This is a proper way to describe both the parent and daughter nuclei from 
the same mean field. Taking {Z, N ) as the even-even nucleus for reference, the parent odd-A nucleus is (Z, N — 1) and 
it decays into the daughter nucleus (Z — 1, N). As an example, to describe the /3 + decay of the odd-neutron parent 
nucleus 73 Kr (Z=36, N=37) into the odd-proton daughter nucleus ,3 Br (Z = 35, A = 38), we use the mean field of 
the even-even nucleus 74 Kr (Z = 36, N = 38). 

In a previous work Jl^] we analyzed the energy surfaces as a function of deformation for all the even-even isotopes 
under study here. For that purpose, we performed constrained HF calculations with a quadrupole constraint |24f| and 
we minimized the HF energy under the constraint of keeping fixed the nuclear deformation. Calculations in this paper 
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are performed for the equilibrium shapes of each nucleus obtained in that way, that is, for the solutions, in general 
deformed, for which we obtained minima in the energy surfaces. Most of these nuclei present oblate and prolate 
equilibrium shapes |l2| that are very close in energy. As we have mentioned, single-particle energies, wave functions, 
and occupation probabilities in the odd nuclei, are obtained from the converged mean fields of the even-even neighbor. 
Therefore, calculations in the odd nuclei are done at the corresponding equilibrium shapes of the even-even generator. 

In this work, besides the odd-nuclei, we also consider a new even-even isotope, 70 Kr, that was not previously 
included in the isotope chain studied in Refs. [i~2|JTT{| . We then show first in Fig. 1 the total energy of 70 Kr as a 
function of the mass quadrupole moment Qq. The results correspond to a constrained HF+BCS calculation with the 
Skyrrne forces SG2 and Sk3. As one can see, we obtain two minima in the energy profile, one is oblate and the other 
is prolate at about the same value of the quadrupole moment. These two minima are separated by about 1 MeV, 
the oblate being the deepest one with the two forces. Thus, 70 Kr follows nicely the trend observed for 72 ’ 74 ’ 76 ’ 78 Kr 
isotopes in Ref. ]Ii|], with an oblate shape favored in a shape coexistent isotope. 


A. Residual interactions 

To describe Gamow-Teller transitions we add to the mean field a spin-isospin residual interaction, which is expected 
to be the most relevant interaction for that purpose. This interaction contains two parts, particle-hole (ph) and 
particle-particle (pp). The ph part is responsible for the position and structure of the GT resonance fdjiq and is 
derived selfconsistently from the same energy density functional (and Skyrrne interaction) as the HF equation, in 
terms of the second derivatives of the energy density functional with respect to the one-body densities |25| . The ph 
residual interaction is finally written in a separable form by averaging the Landau-Migdal resulting force over the 
nuclear volume, as explained in Refs. 00- The coupling strength Xgt is completely determined by the Skyrrne 
parameters, the nuclear radius, and the Fermi momentum. 

The particle-particle part is a neutron-proton pairing force in the = 1 + coupling channel. We introduce this 
interaction in the usual way j^,|l3],^6] 27|, that is, in terms of a separable force with a coupling constant n pp T , which is 
fitted to the phenomenology. Since the peak of the GT resonance is almost insensitive to the pp force, k P q T is usually 
adjusted to reproduce the half-lives J7j] . However, one should be careful with the choice of this coupling constant. 
Since the pp force is introduced independently of the mean field, if k P q T is strong enough it may happen that the QRPA 
collapses, because the condition that the ground state be stable against the corresponding mode is not fulfilled. This 
happens because the pp force, being an attractive force, makes the GT strength to be pushed down to lower energies 
with increasing values of fCg T . A careful search of the optimal strength can certainly be done for each particular case, 
but this is not our purpose in this work. Instead, we have chosen the same coupling constant {k P q T — 0.07 MeV) for 
all nuclei considered here. This value was obtained ||l3|| under the requirements of improving in general the agreement 
with experimental half-lives of even-even isotopes in this mass region while being still far from the values leading to 
the collapse. 


B. Even-even nuclei 


The proton-neutron QRPA phonon operator for GT excitations in even-even nuclei is written as 

= E [ x ™ a t<4 - Y £*«*«*] , 

7TIS 


( 2 . 1 ) 


where a + ( a ) are quasiparticle creation (annihilation) operators, cok are the excitation energies, and , Y^ K the 
forward and backward amplitudes, respectively. From the QRPA equations the forward and backward amplitudes are 
obtained as j2?j 


X UK = 


^K £71 


Y l ° k — 


-1 


UK + €71 


2 X P GT + q„ v M+ K ) - 2 k pp t {£ v M“JL + qX v M%$) 


2x P gt {q^MZ K + q„ v M- K ) + 2 n pp T + qX v M“«_) 


( 2 . 2 ) 

(2.3) 


with e wl/ = E v + E v the two-quasiparticle excitation energies in terms of the quasiparticle energies Ei. M UJK are given 

by 
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where t/s are occupation amplitudes ( u 2 = 1 — t; 2 ) and E-j 
states with spin operators 


+ q*vY£*) , 

(2.4) 

+ q^Y“f) , 

(2.5) 

-ft*) , 

(2.6) 

-0;t) , 

(2.7) 

= Vt/VnY , v k ; , 

(2.8) 


spin matrix elements connecting neutron and proton 


J K 


= {v \(J K \tt) 


(2.9) 


l|. 

by 


Explicit expressions for these matrix elements in the cyl indr ic al h armonic oscillator basis can be found in Ref. 

Note that the expressions for the X and Y amplitudes (2.2)-(2.3) defined here differ from those in Refs. |1 l| . |l2(| 
the contributions due to pp residual interactions that were not include there. 

The solutions of the QRPA equations are obtained by solving first a dispersion relation, which is of fourth order in the 
excitation energies to. Then, for each value of the energy the amplitudes M are determined by using the normalization 
condition of the phonon amplitudes. The technical procedure to solve these QRPA equations is described in detail in 
Ref. 

For even-even nuclei the GT transition amplitudes in the intrinsic frame connecting the QRPA ground state 
|0> (r UK |0> = 0) to one phonon states \ujk) (r+ K |0) = | ujk)), are given by 


(uk\P%\0) = TMf . 


( 2 . 10 ) 


The Ikeda sum rule is always fulfilled in our calculations. 


C. Odd-A nuclei 

The QRPA treatment described in the above subsection is formulated for the excitations of the ground state of an 
even-even nucleus. The GT transition amplitudes connecting the ground state of an even-even nucleus (Oqp state ) 
to all one phonon states in the odd-odd daughter nucleus (2qp states) are given in the intrinsic frame by Eq. ( [2.10 ). 
When the parent nucleus has an odd nucleon, the ground state can be expressed as a lqp state in which the odd 
nucleon occupies the single-particle orbit of lowest energy. Then two types of transitions are possible, which are 
represented schematically in Fig. 2. i) The first type of transitions are phonon excitations in which the odd nucleon 
acts only as a spectator. We call them three quasiparticle transitions (3qp). In the intrinsic frame, the transition 
amplitudes in this case are basically the same as in the even-even case but with the blocked spectator excluded from 
the calculation, 


(/ \Pk \*) 3gp = (“K, lqp \/3 k\ 0, lqp) = tM± k . 


( 2 . 11 ) 


ii) The other type of transitions are those involving the odd nucleon. We call them one quasiparticle transitions (lqp). 
We introduce as usual |||[|27]] phonon correlations to the quasiparticle transitions in first order perturbation, taking 
into account the part of the GT interaction that contains terms linear in the = oYa^ operator. The transition 
amplitudes for the correlated states can be found in Ref. |27|. Here, we give the explicit expression for the transition 
amplitude corresponding to the Kr isotopes, that is, a /3 + -decay of an odd-neutron parent nucleus decaying into an 
odd-proton daughter nucleus. 
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(/ \Pk \ *) i qp — ( 'Kcorr \Pk \ u corr ) — 


C + ^Xgt < C E [( M r ) 2 “«) + (M“ K ) 2 E„(n, cok) 

l ^ K 

-<& E M +" M - K (*, Wif) + (7T, ujf)] 1 

+2k£ t I E uk) - tt, w k )] 


( 2 . 12 ) 


where g's and M's are given in Eqs. ( |2.4|) -(|2.q). In this expression, E n (v,u)K) is in principle given by E„(v,ujk) = 
l/(E 7r — E v —u>k) and a similar expression for E v (rr, u>k) changing v into tt. Nevertheless, in order to avoid accidental 
singularities, we introduce a width d = 0.5 MeV in the same way as it was done for the first time in Ref. || 


E v (v,u>k) = 


E n — E„ — ljk 
{E n -E v - u K ) 2 + d 2 ■ 


(2.13) 


For simplicity we speak of lqp, 2qp, and 3qp states, but it should be clear that we always mean QRPA-correlated 
states. 

Once the intrinsic amplitudes (/ \fd%\ i) are calculated in Eqs. ( 2.11 )-( 2.12 ), the Gamow-Teller strength Bqt for a 
transition — > If can be obtained as 


b gt= e kw iei«>i 2 

= {E W iKil P\ T f K f) \Pp\0Ki) 

P 

+{-l) h ~ Ki {I . _ Ki lp\IfKf) {<t> Kf \f$\4 >kM 

= 8K it Kf [(IiKilO\IfKf) ( 4>K f |/3^| (f>Ki) 

+6 Ki ,i/ 2 (-l) Ii - Ki (h - KillllfKf) (« j> Kf | <t>Ki)] 2 

+^K f ,Ki +1 {IiKill\IfKf) 2 ( <j>K f |/9+i| 4>I<i) 

+5 Kf , Kt -i (IiKd - 1| IfK f ) 2 (<t> Kf |/£ x | 0 Kz ) 2 , (2.14) 


in units of g\/4n. To obtain this expression we have used the initial and final states in the laboratory frame expressed 
in term s of t he intrinsic states \4>k) using the Bohr-Mottelson factorization |^| . 

Eq. (|2.14D can be particularized for even-even parent nuclei. In this case Ii = Ki = 0, If = 1, and Kf = 0,1. 


B gt — 


9a 

4-7T 


+2SK f ,l(<t>K f \/3i\<j>o) }• 


(2.15) 


D. Excitation energies, Qec values, and half-lives 

Concerning the excitation energy of the daughter nuclei to which we refer all the GT strength distributions in 
this paper, we have to distinguish again between the case of even-even and odd-A parents. In the case of even-even 
systems, this excitation energy is simply given by 

E ex [(z,N)^(z-i,N+i)] = u — E n 0 — E„ 0 , (2-16) 

where E %0 and E Uo are the lowest quasiparticle energies for protons and neutrons, respectively. 


5 









In the case of an odd-A nucleus we have to deal with lqp and 3qp transitions. Let us consider here our case of 
an odd-neutron parent decaying by (3 + into an odd-proton daughter. For lqp transitions, since the unpaired neutron 
state is the only neutron state involved in the transition, the excitation energy is 


L^ex.lqp (Z,N — 1)—— 1,N) ~ E n — En 


(2.17) 


On the other hand, in the 3qp case where the unpaired neutron acts as a spectator, the excitation energy with respect 
to the ground state of the daughter nucleus is 


-®'ex,3qp(Z,JV-l)-*(Z-l,JV) — w + inspect — E n 


(2.18) 


This implies that the lowest excitation energy of 3qp type is of the order of twice the neutron pairing gap. In the Kr 
isotopes under study here, A ~ 1.5 — 2 MeV 121. Therefore, all the strength contained in the low excitation energy 
region, below typically 3-4 MeV in the odd-A nuclei studied in this paper must correspond to lqp transitions. Thus, 
the low excitation energy region (below twice the neutron pairing gap) in a /3 + decay from an odd-neutron nucleus, 
basically tells us about the proton spectrum. 

We should also mention that in odd-A nuclei we have considered the transitions to the rotational states as well. 
We consider only those states built on the Kf band heads, which are not forbidden by the selectio n rul es of the 
Gamow-Teller operator in the allowed approximation. This means that we calculate according to Eq. ( 2.14 ) the GT 
strength corresponding to transitions from an initial /, = K \ state to final states characterized by i) If = Kf when 
Kf = Ki + 1, ii) If = Kf and If = Kf + 1 when Kf = Ki, and iii) If = Kf, If = Kp + 1, and If = Kf + 2 when 
K f = Ki- 1. 

Accordingly, we have added to the excitation energies uip a standard Bohr-Mottelson rotational energy (2f| . The 
purpose of the inclusion in our calculation of the transitions to the rotational states is to take into account in a more 
reliable way the fragmentation of the GT strength and the density of excitation energies. Thus, in our calculations we 
estimate the moment of inertia by using a mean value between the two extreme macroscopic models, rigid rotor (rr) 
and irrotational flow (if) model, whose predictions are usually upper and lower boundary values of the experimental 
moments of inertia. 

Qec values are also determined differently depending on the even or odd number of nucleons. The Qec value is 
given by Qec = [M pare nt — A/daughter + m e ] c 2 , where M’s are the nuclear masses (binding energies of the electrons 
have been neglected). For the /3 + -decay of an odd-neutron parent we have 


Qec,(z,n-i)^{z-i,n) — mir — m u + \ 


+ X n ( Z - 2 ,N- 2 ) — A v(Z- 2 .N- 2 ) ~ E, 


v(Z- 2 ,N- 2 ) 


J n(Z- 2 ,N- 2 ) 


E, 


v{Z- 2 ,N- 2 ) j 


(2.19) 


where A is the Fermi level and E is the lowest quasiparticle energy E = y/(e — A) 2 + A 2 . The Q-value for /3 + -decay 
is simply Q 0 + = Q E c ~ 2 m e . 

These expressions can be compared to that corresponding to an even-even nucleus (Z, N ) 


Qec,(z,n)->(z-i,n+i) = — m v + m e + Xtt,(z-2,n) ~ A v ,(z,n) ~ E n f Z -2,N) - E v f Z ,N) ■ 


( 2 . 20 ) 


The half-lives for the /3 + /EC -decay are obtained using the standard definition (see for instance Ref. [^2|), involving 
sums over all possible final states within the Q —window reached in the decay and including standard quenching factors 


The Fermi integrals required for computation of half-lives have been obtained numerically for each value of the 
charge of the daughter nucleus Z and the maximum energy available Wo in /1-decay, as explained in Ref. 301. 


III. RESULTS 

In this Section we present and discuss the results obtained for the GT strength distributions, half-lives, and summed 
strengths in the proton rich Kr isotopes. The results correspond to QRPA calculations with the Skyrme force SG2 
and they have been performed for the nuclear shapes that minimize the HF energy. Before discussing the figures 
we note that the GT strength distributions are plotted versus the excitation energy of the daughter nucleus. The 
distributions of the GT strength in Figs. 3, 5, and 6 have been folded with F = 1 MeV width Gaussians to facilitate 
the comparison among the various calculations, so that the original discrete spectrum is transformed into a continuous 
profile. The distributions in those figures are given in units of g\/I k and one should keep in mind that a quenching of 
the strength is expected on the basis of the observed quenching in charge exchange reactions and spin M 1 transitions 
in stable nuclei, where g s , e fi is also known to be approximately 0.7 g s f lee - Therefore, a reduction factor of about two 
is expected in these strength distributions in order to compare with experiment. This factor is of course taken into 
account when comparison to experimental half-lives, summed strengths, and GT strengths is made. 
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A. Effect of the residual interaction 


We can see in Fig. 3 the effect of the residual interaction treated in QRPA on the uncorrelated calculation (dotted 
lines denoted by MF). We do that on the example of 73,74 Kr for the oblate and prolate shapes that minimize the energy 
obtained with the Skyrme force SG2. The coupling strengths of the ph and pp residual interactions are Xgt = 0.37 
MeV and Kq t = 0.07 MeV, respectively. Both ph and pp residual interactions reduce the GT strength. The residual 
forces produce also a displacement of the GT strength, which is to higher energies in the case of the repulsive ph force 
and to lower energies in the case of the attractive pp force. 

These effects are common to even-even and odd-A isotopes. Nevertheless, we observe that the effect of the residual 
interaction in odd-A nuclei is very small in the low energy tail of the GT strength distribution. This is especially true 
for the pp force. The reason for that can be understood from the fact that the lowest-lying transitions are affected by 
the residual force only through the weak correlations with phonons treated in first order perturbation. 

This feature has also important consequences when one considers the half-lives because they depend only on the 
distribution of the strength below the Qec ~energy window. Therefore, the pp interaction will affect the half-lives 
differently if we deal with even-even or odd-A nuclei. In Fig. 4 we have plotted the half-lives of the 71 - 72 - <3,74j^ r 
isotopes as a function of the strength h P q T . As we can see from this figure the half-lives decrease with increasing 
values of k T q T . but in a different way depending on the isotopes. We can see that the half-lives of odd-A nuclei have 
a very smooth, almost flat, behavior with the pp force, while the half-lives of even-even nuclei present a stronger 
dependence on the pp force. Decreasing of half-lives with h P q T can be understood from the fact that the pp force is 
attractive and therefore tends to concentrate the strength to lower energies. The reason for the smoother decrease of 
half-lives in the odd-A isotopes is related to the smaller effect of the pp force in the low-energy tail of the GT strength 
distribution in odd-A nuclei, previously discussed. 

The optimum value of k P q T to reproduce the half-life depends, among other factors, on the nucleus, shape, and 
Skyrme interaction and a case by case fitting procedure could be carried out. Since the half-lives are practically 
insensitive to this force in odd-A nuclei, this fit could be restricted to even-even nuclei. In Ref. jUI we considered 
this dependence in various even-even isotopes in the mass region A ~ 70 (Ge, Se, Kr, and Sr), and arrived to the 
conclusion that a value of n P q T = 0.07 MeV improves the agreement with experiment in most cases. Thus, we also 
use in this paper the same coupling strength for the pp force. 


B. Comparison of even-even and odd-A GT strengths 

To further illustrate the relative importance of the different types of contributions to the GT strength in odd-A 
nuclei, we show in Fig. 5 the GT strengths (solid) decomposed into their lqp (dotted) and 3qp (dashed) contributions. 
The picture emerging from the analysis of this figure is that the GT strength distributions in odd-A nuclei can be 
divided into two different regions. One is the energy region below twice the pairing gap, where the individual 
excitations are determined by the quasiparticle proton (neutron) energies in the case of an odd-neutron (odd-proton) 
parent nucleus. This region is of relevance for (3 + -decay since it appears within the Q-window. The other region at 
higher energies is dominated by 3qp excitations, where the odd-nucleon acts as a spectator. The strength contained 
here is much larger than in the low energy region because many more configurations are possible but only in the very 
proton rich isotopes this is accessible by (3 + - decay. 

Fig. 6 contains a summary of the results on GT strengths obtained in this work. In this figure we can see our 
HF+BCS+QRPA Gamow-Teller strength distributions predicted by the Skyrme force SG2 for the whole Kr isotopic 
chain including odd-A and even-even nuclei. 

The trend observed in the GT strength distributions of both odd-A and even-even isotopes is similar. We can see 
that the most unstable isotopes have the largest strengths which are also placed at higher energies. As we move into 
the stable isotopes by increasing the number of neutrons, the GT resonance appears at lower excitation energies and 
contains less and less strength. Also the Qec window, represented by the vertical line, becomes smaller and smaller. 

We can observe also the similarity between the even-even and odd-A partners. If we compare the strengths for an 
even-even ( N , Z) nucleus with that of the corresponding odd-A (N — 1, Z) we can see that they are about the same 
once the low energy region of about twice the neutron pairing gap (between 3 and 4 MeV depending on the case) is 
suppressed. 

We can also observe that the different profiles of the GT strength distributions corresponding to the various shapes 
of a given isotope can be used in certain cases as a signature of the nuclear shape. 
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C. Low-energy GT strength and comparison to experiment 


Now we concentrate on the energy region below Qec and discuss the possibilities to discriminate between different 
shapes by /3 + -decay experiments. We also compare with the available experimental data. We perform this detailed 
analysis by plotting not the folded strengths as it was done in the previous figures, but the individual excitations as 
they come from the calculation. The GT strengths have been quenched with the same factor used to calculate the 
half-lives. In Figs. 7-9 we can see these results for the even-even isotopes 72 > 74 > 7 6]£ r) respectively. 

Fig. 7 for ,2 Kr includes the results from our QRPA calculations with the force SG2 with oblate and prolate shapes. 
Below an excitation energy of about 2 MeV we can see that the distribution of the GT strength predicted by the 
oblate or prolate shapes is qualitatively very similar although the prolate shape gives somewhat larger strength. We 
obtain peaks at 0.5 MeV and 1.5 MeV in both cases and they are in agreement with preliminary data on this nucleus 
13 ■ Therefore, it will be hard to distinguish between the two shapes. On the other hand, if we look into the energy 
range from 2 MeV up to Qec , we can see a strong double peak that appears in the oblate case between 2 and 2.5 
MeV. The strength in this region is about three times larger than the strength of the first peak at 0.5 MeV and it 
is almost absent in the prolate case. The appearance or absence of this peak at 2 MeV could be the signature of an 
oblate or prolate shape, respectively. It is also worth mentioning the huge peak appearing in the prolate case very 
close to the Qec value. If it could be seen experimentally, it would be a clear signature in favor of a prolate shape. 
We can see in Table 1 the total measured GT strength below 1.836 MeV |32| compared to our results with the two 
shapes, where oblate shape seems to be favored. 

Fig. 8 shows the results in 74 Kr. Experimental data are from Ref. [[k|. We can distinguish two regions in this 
isotope. At energies below 2 MeV we find that the oblate shape predicts much more strength than the prolate shape 
and the opposite happens beyond 2 MeV. Therefore, measuring the GT strength distribution in 74 I<r up to Qec 
would help to discriminate between the two shapes. If the strength is concentrated below 2 MeV it will correspond to 
an oblate shape, while if it is concentrated between 2 and 3 MeV it will correspond to a prolate shape. A comparison 
of the strength contained below 1 MeV is made in Table 1. The experimental summed strength |l3| lies between the 
predictions of the two shapes but it is closer to the strength produced by the oblate shape. 

In Fig. 9 we show the results for 76 Kr corresponding to the spherical and prolate shapes that minimize the energy 
in this nucleus. Experimental data are from Ref. |m|| . In this case, an isolated single peak in the measured strength 
at very low excitation energy would be the signature of a spherical parent, while a peak close to the Qec limit would 
be the signature of a prolate parent. Comparison with the available data seems to favor the spherical shape. This is 
also true if we compare the total GT strengths contained below 1 MeV, as can be seen in Table 1, where the strength 
generated by the spherical shape is much closer to experiment. 

Next figures 10-13 correspond to the odd-A nuclei 73 75 Kr. In Fig. 10 we can see our results for the oblate and 
prolate shapes of 73 Kr compared to the experimental data from Ref. |2^]. The first thing to mention is that the 
experimental data have been taken differently in the two energy ranges below and above 3.5 MeV. Below 3.5 MeV 
the data have been extracted from direct detection of the gamma rays. Beyond this energy they have been extracted 
from proton delayed detection and have big errors. A new effort to measure experimentally all the energy range up to 
Qec with large efficiency gamma ray detectors is under way at ISOLDE Q. Until these data become available and 
confirm the present measurements, the data measured above 3.5 MeV should be considered as partial because part 
of the strength escapes observation due to the high level density |31|. The summed GT strengths in the two energy 
regions can be seen in Table 1. The experimental data appear between the predictions of the oblate and prolate 
shapes, the oblate one being larger. 

Comparison to the data of the GT strength distributions below 3.5 MeV shows that both oblate and prolate shapes 
produce similar strength in the very low excitation energy range below 0.5 MeV, which is compatible with experiment. 
At higher energies the prolate shape does not generate strength up to 4 MeV while the oblate shape generates two 
bunches at 1 MeV and 3.5 MeV that compare better with experiment. Above 4 MeV, the strength increases and the 
prolate shapes produces two bumps at 4.5 and 6 MeV, while the oblate shape produces a wide peak from 4.5 to 6.5 
MeV with a huge strength at an excitation energy of 5.5 MeV. It is also worth mentioning that the structure of the 
GT transitions in the theoretical calculations can be analyzed taking into account the rotational nature of the final 
states reached by the allowed GT transition. As we have already mentioned, since in this case we are considering the 
transitions K J = 3/2“ —> KJ with Kf = 1/2,3/2, 5/2, when we reach a final state with KJ = 1/2 - , we consider 
the rotational states wi th If = 1/2 - , 3/2 _ , 5/2 _ . The strengths of these states are given by the geometrical Clebsch- 
Gordan factors in Eq. (2.14) and their excitation energies by the rotational energies. Similar arguments apply to the 
transitions to KJ = 3/2 _ , where we consider the rotational states with If = 3/2 _ ,5/2 _ . For KJ = 5/2 _ , the only 
case to account for is If = 5/2 _ . On this basis we can now see that the two peaks below 0.5 MeV in the prolate case 
are the members of a rotational band If = Kf = 3/2“ and If = 5/2~,Kf = 3/2 _ . Similarly we can identify the 
low-lying structure in the oblate case as the members of several rotational bands. This can be seen more easily in Fig. 
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11, where we have plotted the excitation energies in the daughter nucleus 73 Br of the states reached by /3 + - decay. In 
this figure we have also added by dotted lines the first forbidden transitions that correspond to A J = 1 with parity 
change, that is, states l/2 + , 3/2 + , 5/2 + . Although we have not calculated their GT strength, their excitation energies 
have been plotted in the figure. We can see that the agreement with experiment is very reasonable in the oblate case. 
We get concentrations of states at the same energies and the total number of states is similar. In this figure we can 
follow more easily the structure of the bands. For example, in the prolate case we can see that the lowest 3/2“ has a 
5/2“ associated, the next 1/2“ has two states 3/2“ and 5/2“ associated and so on. 

A similar analysis has been done in Figs. 12-13 for 75 Kr. In Fig. 12 we can see the GT strengths in the prolate and 
spherical cases compared to the experimental data from Ref. |Q taken up to 2.5 MeV. The summed strengths can 
be also compared in Table 1. The spherical shape does not produce any strength in the Qec window except a small 
bunch of states around 3 MeV. On the other hand the prolate shape generates strength at the right position although 
a little bit smaller than experiment. This shape also predicts some strength close to Qec, which is not present in the 
spherical case. Fig. 13 shows the experimental low-lying energy spectrum compared to our calculation for the prolate 
shape. In our results we can identify the origin of these excitations. We have a low-lying 3/2 + , 5/2 + , 7/2+ rotational 
triplet and a 5/2 + , 7/2+ doublet. They correspond well with the experimental energies. We should remember that the 
spacing in the rotational energies is determined by the moment of inertia, and we are using a rough estimate. Using 
a moment of inertia a little bit larger, the energy levels become compressed and agree better with the experimental 
spacing. We have also added by dotted lines the energies of first forbidden transitions. 

Experimental information on the decay of 71 Kr is available for the low-lying states J 37 J. In that reference, final 
states for GT decay were reported to the ground state of 71 Br and to excited states at 9 and 207 keV. However, the 
spin-parity assignment of these states, as well as the assignment to the ground state of the parent nucleus 71 Kr are 
still controversial [Q. We can see in Fig. 14 the results of our mean field calculations with the force SG2 for the 
oblate and prolate shapes of 72 Kr. The vertical axis is the neutron occupation probability and the horizontal one 
is the neutron single-particle energy. The deformed single-particle states are labelled by A' 71 ’. We can see that the 
picture is compatible with several spin and parity assignments for the odd-neutron isotope 71 Kr. If we consider the 
neutron level closest to Fermi level in Fig. 14, the spin-parity assignment would be 9/2 + for the oblate shape and 
3/2“ for the prolate shape. On the other hand, the assumption taken in Ref. [^7) was 5/2 _ , which is also close to 
Fermi level. We therefore consider these two possibilities for each shape. 

We show in Fig. 15 the results for the GT strengths from our calculations assuming the above mentioned possibilities 
for the spin and parity of the parent nucleus, 5/2~ as it was taken in Ref. |57|, as well as 9/2 + in the oblate case and 
3/2“ in the prolate case. We can also see the spin and parity of each GT excitation. 


D. Qec values and half-lives 

Experimental Qec and Ti/ 2 values are plotted as circles in Fig. 16 for the Krypton isotopes considered in this work. 
The data have been taken from Ref. |2l]], except for the isotopes 70,71 Kr, where we have used the more recent data 
from Ref. j37j. These values are compared to our theoretical results represented by the vertical lines. The extreme 
values of these vertical lines correspond to the results obtained from the various shapes using the forces SG2 and Sk3. 
This has been done in order to have a better idea of the theoretical spread of the results. The agreement is in general 
good for both Qec and half-lives. This is notorious because there is a very large range of variation (seven orders of 
magnitude) for half-lives. Qec values are well reproduced with the exception of ,0 Kr, where we obtain a value below 
experiment. 

We can also see that the calculations fail to account for the half-lives of the most unstable nuclei 70 Kr and 71 Kr. 
Nevertheless, it should be mention that the calculations correspond to GT transitions neglecting possible contributions 
from Fermi transitions. But Fermi transitions might play a significant role to obtain the total decay rates in nuclei 
where N ~ Z. Indeed, we have calculated the Fermi strength distributions and evaluated the corresponding half-lives. 
The calculation of Fermi transitions follows closely that of Garnow-Teller transitions. One should simply use an 
isospin-isospin residual interaction with a coupling strength derived from the Skyrrne force as shown in Ref. flT| and 
replace the spin matrix elements between proton and neutron states given in Eq. by the overlaps < v\t: >. Of 

course, only K = 0 components will survive. 

We find that the half-lives corresponding to Fermi transitions are negligible as compared to the GT half-lives except 
in the isotopes 70 > 71 ,72,73 w j lere they are comparable and therefore one should not neglect their contribution. We 
can see in Fig. 16, the new results for the total half-lives once the contribution from Fermi transitions have been 
included. They are given by the vertical lines located at the right side in the half-lives of the ' ( V 71 - 72 ’ 73 Kr isotopes. 
The agreement with experiment improves in the cases of 70,71 Kr but we are still above the experiment. Nevertheless, 
it should be mentioned that if we use the experimental Qec value in '°Kr, instead of the Qec result from our 
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calculations, we get agreement with the experimental half-life. Thus, the discrepancy found in the half-life of 70 Kr is 
entirely due to the discrepancy in the Qec value. 

We can observe in Fig. 16 a nice trend that can be followed by the auxiliary dotted lines joining the experimental 
data of even-even isotopes on one hand and the odd-A on the other. The growing Qec trends in the two lines of Fig. 
16 are dictated by the increase of (A w — A^) as one approaches the pro ton d ri p line . This behavior can be qualitatively 
understood in the case of Qec from their expressions given by Eqs. ( |2.19| )-( 2.20 ). As we have mentioned earlier, the 
only exception to this trend appears in the experimental Qec value of 70 Kr, which is clearly above the theoretical 
value. Since in this case, the daughter nucleus 70 Br is odd-odd N = Z, one may argue that the observed jump from 
the regular trend in Qec and Tj / 2 is due to an extra neutron-proton binding, which is not taken into account in the 
present calculation. 


IV. CONCLUSIONS AND FINAL REMARKS 

We have applied a selfconsistent deformed HF+BCS+QRPA formalism with density-dependent effective Skyrme 
interactions to the description of the /3-decay properties of proton-rich odd and even Kr isotopes. This approximation 
has the appealing feature of treating the excitations and the ground state in a selfconsistent framework with basically 
no free parameters. Our spin-isospin residual interaction contains a particle-hole part, which is derived selfconsistently 
from the Skyrme force, and a particle-particle part, which is a separable force representing a neutron-proton pairing 
force. 

We have analyzed the similarities and differences in the treatment and in the results of even-even and odd-A nuclei. 
The low-lying GT response of odd-A nuclei is generated by transitions involving the state of the odd nucleon (lqp 
transitions). Thus, the excitation spectrum up to twice the pairing gap parameter of neutrons (protons) in the /3 + - 
decay of and odd neutron (proton) parent, gives information on the proton (neutron) states. When this low energy 
strength in the odd-A nuclei is removed, the resulting GT strength distribution is very similar to that of the even-even 
neighbor displaced to higher excitation energy ( E ex ~ 2A). 

In odd-A deformed nuclei, for each allowed intrinsic GT transition, one has a set of transitions to rotational states 
with decreasing strength at higher energies. This is a characteristic of odd-A deformed nuclei, which is not present in 
the spherical limit. The GT strength corresponding to transitions to states in a given band is reduced by geometrical 
factors involving angular momenta. The energy separation depends on the angular momentum of the transition with 
a global scale determined by the moment of inertia. 

We have found a reasonable agreement with available experimental data. We have also discussed what can be 
learned from future comparison of our results with experimental data from ISOLDE that are expected very soon. 
This comparison would be interesting for several reasons: i) From our study of the dependence on the shape of the 
GT strength distributions we conclude that information on the shape of the parent nucleus can be gained when data 
in the whole Qec window become available. We have identified particular narrow regions in the excitation spectra of 
various nuclei where data could be more conclusive on the nuclear shape; ii) Comparison to data would be desirable 
before considering further theoretical refinements, such as the effects of the continuum, the neutron-proton pairing 
at the mean field level or extensions of the RPA; iii) An important point will be to check whether our calculations 
produce the same level of agreement in all the isotopes or some special characteristics can be found in N = Z nuclei. 
For the moment, an interesting feature found is the deviation between theory and experiment in the half-life of '°Kr. 
From the present calculation, this deviation can be attributed to the difference between experimental and theoretical 
Qec values, which in turn can be a signature of an extra binding in the N = Z odd-odd daughter nucleus. 
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Table 1. Com par ison of the GT strengths contained below some given excitation energy between experimental 
measurements ( 32 for ' 2 Kr, HH for ,3 Kr, for ,4 Kr. J3^] for 75 Kr, and Q] for 76 Kr), and theoretical calculations. 



exp 

oblate 

prolate 

72 Kr ( E ex < 1.836 MeV) 

0.5 ± 0.1 

0.5 

0.8 

73 Kr ( E ex < 3.5 MeV) 

0.10 ± 0.02 

0.19 

0.06 

73 Kr (4.0 MeV < E ex < 6.5 MeV) 

0.83 ± 0.60 

0.98 

0.42 

74 Kr (E ex < 1 MeV) 

0.20 

0.23 

0.12 

76 Kr ( E ex < 2.2 MeV) 

0.08 

0.00 (sph) 

0.05 

76 Kr ( E ex < 1 MeV) 

0.16 

0.18 (sph) 

0.05 
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FIG. 1. Total energy of 70 Kr as a function of the mass quadrupole moment Qo obtained from a constraint HF+BCS 
calculation with the Skyrme forces SG2 (solid) and Sk3 (dashed). The distance between ticks in the vertical axis corresponds 
to 1 MeV but the origin is different for the two forces. 
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FIG. 3. Gamow-Teller strength distributions \g\l^\ in ,3 ’ ,4 Kr isotopes plotted versus the excitation energy of the cor¬ 
responding daughter nucleus 73 ' ,4 Br, respectively. The calculations are performed in HF+BCS approximation (dotted lines 
denoted by MF) and in QRPA including only the ph residual interaction (dashed) and including both ph and pp residual 
interactions (solid). 
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FIG. 5. Decomposition of the total GT strength distribution into their lqp and 3qp contributions (see text). 
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FIG. 6. Gamow-Teller strength distributions [g\/kn\ as a function of the excitation energy of the daughter nucleus [MeV]. 
The results correspond to the Skyrme force SG2 in QRPA for the various shapes of the even-even and odd-A Krypton isotopes. 
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FIG. 11. Experimental and calculated decay schemes for 73 Kr. 
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FIG. 13. Same as in Fig. 11 for the decay of 75 Kr. 
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FIG. 14. Neutron single-particle energies and occupation probabilities of the intrinsic states K n calculated with the force 
SG2 for the two shapes, oblate and prolate, that minimize the energy in 72 Kr. A„ are the neutron Fermi energies. 
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FIG. 15. Energy distribution of the Gamow-Teller strength in 71 Kr assuming different spins and parities of the ground state 


in the parent nucleus. 


28 




70 71 72 73 74 75 76 77 

Kr 


FIG. 16. Experimental Qbc values (top) and half-lives (bottom) for the Kr isotopes are given by circles. The solid vertical 
lines correspond to our theoretical QRPA calculations. The lengths of the vertical lines indicate the different results we obtain 
from using the forces SG2 and Sk3 and the possible shapes in each isotope. 
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